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Abstract 

The RNA chaperone Hfq acts as a central player in post-transcriptional gene regulation in several Gram-negative Bacteria, 
whereas comparatively little is known about its role in Gram-positive Bacteria. Here, we studied the function of Hfq in 
Bacillus subtilis, and show that it confers a survival advantage. A comparative transcriptome analysis revealed mRNAs with a 
differential abundance that are governed by the ResD-ResE system required for aerobic and anaerobic respiration. 
Expression of resD was found to be up-regulated in the hfq~ strain. Furthermore, several genes of the GerE and ComK 
regulons were de-regulated in the hfq~ background. Surprisingly, only six out of >100 known and predicted small RNAs 
(sRNAs) showed altered abundance in the absence of Hfq. Moreover, Hfq positively affected the transcript abundance of 
genes encoding type I toxin-antitoxin systems. Taken the moderate effect on sRNA levels and mRNAs together, it seems 
rather unlikely that Hfq plays a central role in RNA transactions in Bacillus subtilis. 
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Introduction 

In response to environmental cues, Bacteria synthesize small 
regulating RNAs (sRNAs). Many sRNAs act at the post- 
transcriptional level through base-pairing with target mRNAs 
[1]. In Escherichia coli, the sRNA-mRNA interactions often require 
the RNA chaperone Hfq. Hfq belongs to the Sm and Sm-like 
family of proteins [2,3]. It interacts with sRNAs and mRNAs, and 
assists in sRNA-mRNA annealing. Hfq-mediated riboregialation 
can result in post-transcriptional gene activation or silencing [3]. 
In addition, binding of Hfq to sRNAs can affect their stability [4- 
7]. Comparative genome analyses revealed that approximately 
half of the sequenced bacterial genomes contain at least one gene, 
encoding a Hfq orthologue [2,8-9]. Throughout bacterial species, 
Hfq proteins comprise a conserved core but differ in length of their 
C-terminal extensions, with the y- and fi-Proteobacteria displaying 
the longest C-termini [9] . In contrast, the Hfq orthologs of Gram- 
positive Bacteria lack extended C-termini [2,8]. A c-terminaUy 
truncated E. coli Hfq variant comprising the conserved core (aa 1 - 
65) was unable to support translational autoregulation of the hfq 
gene as well as translational activation of the rpoS gene by the 
sRNA DsrA at low temperature [10,11]. However, Hfq variants 
containing N-terminal 69 and 72 aa, respectively, were proficient 
in mediating riboregulation of other genes [1 2] . In light of these 
results it seems possible that some Hfq-mediated regulatory events 
require the extended C-terminus, whereas it seems to be 
dispensable for others. 



The Bacillus subtilis hfq gene (first annotated asjmaH) encodes a 
73 amino acid protein (Hfq^,), forms homohexamers and was 
shown to bind to a RNA aptamer [13]. Hfq^, was unable to 
functionally replace the E. coli protein in vivo [10]. Similarly, the 
presence of Staphylococcus aureus Hfq did not rescue Hfq-related 
phenotypes of a Salmonella typhimurium hfq mutant [14]. Despite of 
producing Hfq.sa at levels comparable to Hfq^, in wild-type 
Salmonella, the complemented strain behaved as a hfq null mutant 
in three sRNA-mediated regulatory responses. In addition, in 
contrast to E. coli Hfq, Hfq^, did not stimulate annealing of 
complementary oligonucleotides, suggesting that the evolutionarUy 
conserved core may not be sufficient to support base-pairing of 
two RNAs [11]. 

The B. subtilis SRI sRNA acts by base-pairing and inhibits 
translation of ahrC mRNA, encoding a positive regulator of the 
arginine catabolic operons rocABC and rocDEF [15]. Although B. 
subtilis Hfq was shown to bind to both, SRI sRNA and ahrC 
mRNA in vitro, it did not stabilize the SRI sRNA in vivo, and was 
not required for duplex formation between the two RNAs. 
However, the expression of a translational ahrC reporter gene 
fusion was reduced in the absence of Hfq, indicating that the 
protein is required for efficient expression of ahrC mRNA in vivo 
[15,16]. So far, no evidence for an involvement of Hfq in sRNA- 
mediated regulation has been obtained in B. subtilis [15-20], 
Nevertheless, a deep sequencing approach of Hfq-associated 
RNAs revealed a subset of sRNAs, antisense RNAs, antitoxin 
transcripts, and many mRNA leaders that associate with Hfq [21]. 
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In Gram-positive Bacteria, a phenotype for Hfq has only been 
described in S. aureus and Listeria monocytogenes. Liu et al. [22] 
showed that the deletion of hfq caused increased carotenoid 
pigmentation and a decreased toxicity of S. aureus, the latter of 
which might be attributed to the differential expression of more 
than 100 genes, most of which are related to stress response and 
virulence of S. aureus. However, for the RNAIII/ spa model system 
it was shown that Hfq did not affect sRNA-mRNA annealing 
although it specifically bound to RNAIII and spa mRNA in vitro 
[23]. Furthermore, Hfq did not affect RNAIII or spa mRNA 
quantities in vivo [24] . In addition, Hfq is apparently not required 
for other targets of RNAIII, sal 000 and coa mRNAs [24,25]. In 
accordance, Hfq had no major impact on the levels of 1 1 sRNAs, 
characterized in S. aureus [26]. In E. coli, Hfq-dependent sRNA- 
mRNA duplexes are formed between relatively short RNA base- 
pairing sequences, whereas sRNA-mRNA duplexes in S. aureus are 
generally more stable [27], which could explain the dispensability 
for Hfq-like RNA chaperones in this organism. 

In L. monocytogenes Hfq was shown to be important for tolerance 
to osmotic and ethanol stress, for long-term survival under amino 
acid-limiting conditions and for pathogenicit)' in mice [28]. Co- 
immunoprecipitation of Hfq-associated RNAs identified three 
novel sRNAs, LhrA-C [29] . .Although Hfq interacted with all three 
sRNAs in vitro, only LhrA sRNA showed a reduced stability in a Z,. 
monocytogenes hfq~ strain [29]. Further studies revealed that LhrA 
most likely acts by regulating the translation of three mRNAs, 
lmo0850, lmo0302 and chiA through an anti-sense mechanism that 
depends on Hfq [30,31]. 

In Streptococcus pneumoniae, close to hundred sRNAs have been 
identified [32-37]. However, most of their target genes remain 
unidentified. Some of the S. pneumoniae sRNAs are part of the 
regulon governed by the two-component system CiaRH [32]. 
Computational target predictions and verification by reporter gene 
assays identified six genes that are controlled by five redundant 
csRNAs (da-dependent .small RNAs): sprOOSl, spr0371, spr0159, 
spr0551, sprl097 and spr2043 [37]. Thus, sRNA-mediated post- 
transcriptional regulation seems to occur in S. pneumoniae, although 
the family of Streptococcaceae does not contain a hfq homologue 
[9]. 

In this study we have performed comparative transcriptional 
profiling of a B. subtilis wild-type strain and an isogenic hfq deletion 
mutant using an RNAseq approach. In contrast to studies in 
Gram-negative bacteria [38,39] the absence of Hfq had no global 
effects on the transcriptome. Hfq rather impacted on distinct 
regulons including the (i) ResD-ResE signal transduction .system 
required for aerobic and anaerobic respiration [40] , (ii) the GerE 
regulon encoding sigma-K-dependent late spore coat genes [41], 
as well as (iii) the ComK regulon, required for the regulation of 
competence and DNA uptake [42]. Surprisingly, only six known 
and predicted B. subtiiis sRNAs showed an altered abundance. 
Furthermore, 3 of the down-regulated RNAs in the hfq~ strain 
belong to the class I toxin-antitoxin (TA) systems and represent the 
toxin encoding mRNAs. In the absence of Hfq, the levels of their 
convergentiy transcribed antitoxin/ antiscnse RNAs were also 
reduced to a similar extend as the toxin mRNAs. Thus, Hfq 
seems to affect both the toxin and antitoxin transcripts in a positive 
manner. 

Results and Discussion 

The hfq Gene Confers a Survival Advantage to 6. subtilis 
in Stationary Phase 

To monitor ffq expression throughout growth, B. subtilis 168 was 
cultivated in minimal CS-glucose medium. The hfq mRNA levels 



were determined by primer extension analysis at different times 
during growth (Figure lA). In addition, Hfq synthesis was assessed 
in strain GP1067, encoding a hfq::flag fusion gene, transcription of 
which is driven by the authentic promoter. In minimal CS-glucose 
medium hfq was expressed in B. subtilis 168 during early 
logarithmic growth and expression only slightly incrc-ascd in late 
stationary phase (Figure lA). Similarly, no gross difference in the 
Hfq levels was observed in strain GP1067 throughout growth 
(Figure lA). The latter result is at variance with a recent study 
[21], reporting that the synthesis of chromosomally encoded Hfq- 
Flag increasc-d upon transition to stationary phase. 

To test whether Hfq impacts on the viability of & subtilis, strain 
168A^ was co-cultivated with the parental wild-type strain and 
strain GP1067, respectively, and growth was monitored in CS- 
glucose minimal medium at 37°C over 10 days. As shown in 
Figure IB, the CFU of the \&Q\kfq strain started to decline at day 
two after inoculation. On day three, the CFU of strain l68Ahfq 
was reduced to ~ 1 % when compared with the wild-type strain or 
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Figure 1. Expression of the hfq gene and requirement of Hfq 
for survival of B. subtilis. (A) Bacillus subtilis strain 1 68 wt and GP1 067 
were grown in CS-glucose medium at 37 C. The hfq mRNA and 5S rRNA 
(internal control) levels were assessed by primer extension analyses in 6. 
subtilis strain 168 (upper two panels). Strain GP1067 was used to 
monitor Hfq and ribosomal S2 (internal control) levels by quantitative 
western-blotting (lower two panels). Samples for RNA extraction and for 
western-blot analyses were withdrawn throughout growth at the OD^oo 
values indicated on top. Only the relevant sections of the autoradio- 
graphs and immunoblots are shown. (B) The 6. subtilis strains 168 wt 
(triangles) and ISSAhfq (squares) were co-cultivated in CS-glucose 
medium. Likewise, 6. subtilis strains GP1067 (diamonds) and ^68Ahfq 
(circles) were co-cultivated in CS-glucose medium. Growth was 
monitored over 10 days by scoring the CFU. Immunodetection of 
Hfq-Flag (lower panel) in strain GP1067 was performed at the days 
indicated. Only the relevant section of the immunoblot is shown. 
dol:1 0.1 371/journal.pone.0098661 .gOOl 
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with strain GP1067. As shown for other Bacteria [10,28,43], this 
experiment indicated that Hfq might play a role in survival of B. 
subtilis cells during stationary phase and prompted us to compare 
the transcriptomes of B. .subtilis strains 168 and 168AA/^ by 
RNAseq during logarithmic growth and at the onset of stationary 
phase. 

Transcriptional Profiling of 6. subtilis ^68Ahfq 

Total RNA was isolated from cultures of B. subtilis 1 68 vvt and 
168A^ grown in CS-glucose minimal medium at 37°C to an 
ODgoo of 0.7 and 2.0, respectively. The growth rate was 
indistinguishable until an ODgoo of ~0.7, whereas the growth 
rate of the hfq strain was somewhat reduced at higher ODgoo 
values (Figure SI). All genes, annotated in the NCBI database 
were included in the differential gene expression analysis. A p- 
value (adjusted for multiple testing) of 0. 1 was set as threshold for 
significance. 

When compared with the wild-type strain 168, 68 and 88 
mRNA transcripts showed a differential abundance in strain 
168Ahfq during logarithmic growth and in early stationary phase, 
respectively (Figure 2A and B). Among them, the level of 43 
mRNAs was altered during both growth phases. In logarithmic 
growth phase 30 mRNAs were down-rc-gulated and 38 were up- 
regulated (Table SI in File SI), whereas in early stationary phase 
22 were down- and 66 were up-regulated (Table S2 in File SI). A 
functional classification of these genes revealed that a significant 
number of transcripts belong to the ResD, Rex, GerE and ComK 
regulons (Figure 3A and B). 

Expression of Genes Involved in Anaerobic Respiration is 
Enhanced in the Absence of Hfq 

Notably, 16.2'!^) and 28.4'X) of the mRNAs with significantiy 
altered abundance- in logarithmic and early stationary phase, 
respectively, are governed by the ResD-ResE two component 
system required for aerobic and anaerobic respiration [40]. The 
abundance of transcripts of the resABCDE operon, as well as the 
abundance of transcripts of many ResD-regulated genes was 
increased in the hfq strain (Tables SI and S2 in File SI). ResE is a 
sensor kinase that, upon autophosphorylation, phosphorylates the 
response regulator, ResD. Subsequendy, phosphorylated ResD 
activates transcription at target promoters. The resD and resE genes 
constitute an operon with the three upstream genes, resABC (Figure 
S2). The resD/E genes are transcribed from a resDE-fj^tciSx, 
promoter and the resA operon promoter, the latter of which is 
dependent on ResD/E. The ResD/E system governs anaerobic 
respiration by regulating the anaerobic regulator fnr and the 
nasDEF genes, which encode subunits of nitrite reductase [44—47] 
(Figure S2). The genes Jnr, nasD and nasE displayed an increased 
transcript abundance in the hfq~ strain (Tables SI and S2 in File 
SI). 

The genes cydABCD, IctP-ldh andywcj, which are regulated by 
the transcriptional regulator Rex, were as well up-regulated in the 
absence of Hfq (Tables SI and S2 in File SI). The Rex repressor 
controls expression of the respiratory oxidase cytochrome bd 
{cydABCD), the NADH-Hnked fermentative lactate dehydrogenase 
{IctP-ldh), a type II NADH dehydrogenase (yjlC-ndh) and a putative 
nitrite transporter (ywcj) [48-50]. The redox-sensing Rex protein 
indirectly responds to changes in oxygen availability and 
transcription of the Rex-repressed genes is activated when oxygen 
is limiting [48] . In the absence of oxygen, ResD-dependent genes 
are activated as well. As the transcript abundance of several genes 
of the ResD and Rex regulons was increased in the hfq~ mutant, it 
seems likely that Hfq influences the adaptation of B. subtilis to 
anaerobic growth. This could be mediated by the ResD-ResE 



system as it occupies an early stage in the regulatory pathway 
governing anaerobic respiration. 

To verify up-regulation of the resD/E transcripts in the absence 
of Hfq the B. subtilis strain 168 and the isogenic strain 168AA/^ 
were transformed with plasmids encoding transcriptional resA-gfp 
and resD-gfp - and translational resA::gfp and resD::g/p fusion genes, 
respectively. The strains were grown in CS-glucose medium at 
37°C to an ODgoo of 0.7 when fluorescence of the cells was 
measured. For both genes, resA and resD, transcription was 
enhanced in the hfq background when compared with the hfq'^ 
strain (Figure 4A and B), which was also mirrored by the increased 
translation of the resAr.gfp and resD::gfp fusion genes (Figure 4A and 
B). As up-regulation of resD transcription results in elevated levels 
of resA [47], we assume that Hfq negatively impacts primarily on 
the transcript levels of resD. 

In E. coli, Hfq acts as a RNA chaperone, which mediates sRNA- 
mRNA interactions and thereby regulates the expression of target 
mRNAs [3]. As Hfq^, affected the transcript abundance of the 
resABCDE operon in a negative manner (Figure 4A and B), we 
hypothesized that the transcripts of the resABCDE operon could 
represent targets that are post-transcriptionally regulated by Hfq 
and sRNAs. A bioinformatic approach was used to identify sRNAs 
that could possibly interact with the transcripts of the resABCDE 
operon. Verified and predicted sRNAs oi Bacillus subtilis were used 
as a source for the analysis [51,52]. Three sRNA candidates 
(sRNA 13 [51], FsrA [19], SurA [53], Figure 5A) with the 
potential to base-pair adjacent to the start codon of resA were 
identified. However, none of the sRNAs showed significandy 
altered levels in the RNAseq analysis, and we did not observe a 
difference in abundance of these sRNAs in the absence of Hfq in 
Northern-blot analyses (Figure 5B). The unchanged levels of FsrA 
observed in the hfq mutant (Figure 5B) are at variance with recent 
data by Dambach et al. [21], who observed reduced levels of FsrA 
in the absence of Hfq. Transcription of FsrA is repressed by the 
ferric uptake regulator protein Fur [17]. Interestingly, resA mRNA 
levels were found to be decreased in a B. subtilis fur mutant that has 
elevated levels of FsrA. However, the levels oiresA were restored to 
wild-type levels in fur/fbpAB or fur/fbpABC mutant strains, 
suggesting that the potential regulation of resA by FsrA is mediated 
by the FbpA, B and C proteins rather than by Hfq [17]. FpbA,B 
and C are small basic proteins, which have been proposed to 
function as dedicated Fur-regulated RNA chaperones [17]. 

Hfq Affects the GerE and ComK regulons 

17,6% and 10,2% of all mRNAs with significantiy altered 
abundance in logarithmic and early stationary phase, respectively, 
belonged to the GerE regukm (Tables SI and S2 in File SI). The 
regulon is controlled by the transcriptional regulator GerE and 
encodes a subset of sigma-K-dependent late spore coat genes [41]. 
The transcript levels oi gerE and many genes governed by GerE 
are decreased in the hfq~ strain (Tables SI and S2 in File SI). In 
opposite, many genes, which were up-regulated in the presence of 
Hfq belonged to the ComK regulon, required for the regulation of 
cell competence and DNA uptake (Tables SI and S2 in File SI) 
[42,54]. Prolonged nutritional stress can result in the development 
of competence and sporulation [54,55]. Towards sporulation, the 
individual cells can opt for the diflFerentiated state of competence, 
triggered by the competence master regulator ComK [54-56]. In 
this state the cell can take up exogenous DNA from lysed cells. As 
Hfq impacted inversely on the ComK and GerE regulons, it could 
be involved in regulating steps towards sporulation. 
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Figure 2. DEseq analysis of B. subtilis 168 wt versus leSA/j/t/. (A) Differential abundance of transcripts in strain 168 wt and in strain 168A/)/q 
during logarithmic growth (OD6oo = 0.7) (B) Differential abundance of transcripts in strain 168 wt and in strain 168Ah/q in early stationary phase 
(OD5oo = 2.0). Each dot represents one transcript. The log2 fold-change is plotted against the mean expression level for each transcript. Red dots 
represent transcripts whose abundance is significantly altered (p-value adjusted for multiple testing <0.1). 
doi:1 0.1 371 /journal.pone.0098661 .g002 



sRNA Abundance is Hardly Affected by Hfq 

In E. coli, Hfq plays a central role in sRNA-mediated post- 
transcriptional regulation of target mRNAs, and was shown to 
affect the stability and half-life of several sRNAs [3] . Therefore, it 
was somewhat surprising that only a few B. subtilis sRNAs showed 
a significant change in abundance upon deletion of the hfg gene. 
Only six out of the > 1 00 predicted and confirmed sRNAs of B. 
subtilis [51,52] were found to be diflFerentially abundant (Tables SI 
and S2 in File SI). All other sRNAs did not display a significant 
change in abundance in the absence of Hfq as the corresponding 
p-values (adjusted for multiple testing) were above the threshold set 
for significance (>0.1). However, five of the sRNAs exhibited a 
decreased abundance in the absence of Hfq, whereas the level of 
one sRNA (sRNA 471) was increased (Tables SI and S2 in File 

Similarly to this study, only 22 B. subtilis sRNAs were identified 
by Co-IP with Hfq^,, specific antibodies [2 1] . Correspondingly, a 
Co-IP approach in Listeria monocytogenes revealed only three sRNAs 
[29], whereas more than 60 sRNAs have been predicted or 
confirmed in L. monocytogenes [57-59]. In contrast, RNA Co-IP 
experiments with epitope-tagged Hfq in Salmonella resulted in the 
detection of about half of the sRNA complement in this organism 



[60]. In contrast to Enterobacteriaceae it seems therefore rather 
questionable whether Hfq plays a central role in sRNA-mediated 
post-transcriptional regulation in Firmicutes. 

Hfq Affects the Abundance of Toxin/Antitoxin RNAs 

Among the mRNAs that exhibited a signilicantiy reduced 
abundance, 3 mRNAs encode the toxin genes bsrE, bsrG and bsrH, 
which belong to type I toxin-antitoxin (TA) systems. In type I TA 
systems the antitoxin is an antisense RNA that is convergentiy 
transcribed with the toxin encoding mRNA [61]. Interestingly, the 
antisense RNAs of the bsrE, bsrG and and txpA TA systems were 
also enriched in Co-IP studies with Hfq^j [21]. To verify the 
reduced levels of the three toxin mRNAs in the absence of Hfq, 
strain 168 and the isogenic \68 Ahfq strain were grown in CS- 
glucose medium at 37°C and samples for RNA extraction were 
withdrawn at an ODgoo of 0.7 and 2.0. As verified by Northern- 
blot analyses, the levels of the hsrE, bsrG and bsrH toxin mRNAs as 
well as their corresponding antitoxin asRNAs as-BsrE, SR4 and 
as-BsrH were reduced in the absence of Hfq (Figure S3). In line 
with our study, the transcript levels of the type I TA systems bsrE/ 
asBsrE, AjrG/ SR4 and the RatA antitoxin were reported to be 
decreased in a. hfq strain [21]. In addition, our study revealed a 
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Figure 3. Functional classification of mRNAs with altered abundance in B. subtilis strain "iSS^hfq. (A) Gene transcripts significantly altered 
during logarithmic phase. (B) Gene transcripts significantly altered during early stationary phase. Genes are classified either according to their regulon 
affiliation (for the regulons ResD, Rex, GerE and ComK) or according to gene ontology (restricted to the terms within the domain "biological 
process"). 

doi:1 0.1 371 /journal.pone.0098661 .g003 



decrease in the abundance of the transcript levels of bsrH/ asBsrH 
in the absence of Hfq (Figure S3). However, Hfq neither affected 
the half-life of the A.trG/SR4 and txpA/RatA RNAs nor was it 
required for the function of the TA systems [18,19]. Thus, it is 
tempting to speculate that Hfq positively affects transcription of 
these TA transcripts rather than their stability. 

Concluding Remarks 

The comparative transcriptome analysis revealed mRNAs with 
a dififerential abundance, which are governed by the ResD-ResE 
system as well as several genes of the GerE and ComK regulons. 
Hence, it will be interesting to focus in further studies on the role 
of Hfq in adaptation to anaerobic growth and whether it affects 
competence and sporulation processes. The mRNAs/asRNAs of 
several type I TA systems were reduced in the absence of Hfq. As 
these TA systems have been hypothesized to confer an advantage 
to cells during stress conditions, it seems further interesting to ask 
whether there is a link between the reduced viability of the hfq— 
strain (Fig. IB) and the reduced expression of the type I TA 
systems. 

The comparison of the transcripts displaying an altered 
abundance in the absence of Hfq (Tables SI and S2 in File SI) 
with mRNAs that associated with Hfq^, [21] revealed significant 



differences. Approximately 150 Hfq-bound mRNAs fragments 
were detected in the Co-IP screen. However, only five mRNAs 
were identified in both approaches. Notably, of the five genes 
identified in both studies, two {ctaD and sboA) are regulated by the 
ResD-ResE signal transduction system, while the other three genes 
(yxiE, yebD and jdbL) are of unknown function. It is conceivable 
that the altered levels of many transcripts in the absence of Hfq 
result from a perturbation of the expression levels of a limited 
number of regulatory genes, and thus from indirect effects. In 
agreement with this is the observation that the abundance of 
transcripts of distinct regulons was affected (Fig. 3 A and B). 

Materials and Methods 

Bacterial Strains, Growth Conditions and Transformation 

The strains used in this study are listed in Table S3 in File SI. 
All B. subtilis strains used in this study are derivatives of the wild- 
type strain 168. Escherichia coli TOP 10 (Invitrogen) was used for the 
cloning experiments. B. subtilis was grown in Luria Bertani (LB) 
medium or C-minimal medium with succinate [62] in the presence 
of Q.y/a glucose (CS-glucose medium) supplemented with chlor- 
amphenicol (5 |ag/ml), spectinomycin (150 (ig/ml) or phleomycin 
(2.5 |a.g/ml), where appropriate. 
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Figure 4. Down-regulation of A'e5>4/I7 expression in tlie presence of Hfq. (A) Fluorescence conferred by plasmid pADresAts (1) borne 
transcriptional resA-gfp fusion and by the plasmid pADresAtI (2) borne translational resAr.gfp fusion, respectively, was determined in S. subtilis strains 
168 wt (black bars) and 168Ah/q (gray bars). (B) Fluorescence conferred by the plasmid pADresDts (1) borne transcriptional resD-gfp fusion and by the 
plasmid pADresDtl (2) borne translational resD::gfp fusion, respectively, was determined in B. subtilis strains 168 wt (black bars) and 168A/iftj (gray 
bars). Error bars represent standard deviations. 
doi:10.1371/journal.pone.0098661.g004 
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Figure 5. TKie abundance of putative sRNA regulators of resA is independent of Hfq. (A) Partial sequence of the 5'-untranslated region of 
resA and of the initial coding region. The —10 and —35 regions of the resA operon promoter are highlighted in blue. Arrow: start site and direction of 
transcription. The start codon (ATG) of resA is indicated in red. Predicted base-pairing interactions of sequences of FsrA, SurA and sRNA 13 with the 
rbs of resA mRNA. (B) The levels of FsrA, SurA, sRNA 13 and 5S rRNA (loading control) were determined by Northern-blot analyses as described in 
Materials and Methods. Bacillus subtilis strains 168 wt (wt) and ISSAhtq (A) were grown in CS-glucose medium at 37 'C. Samples for RNA extraction 
were withdrawn at the ODgoo values indicated on top. Only the relevant sections of the autoradiographs are shown. 
doi:1 0.1 371/journal.pone.0098661 .g005 
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The B. subtilis strain l68A/ifq was constructed by transformation 
of chromosomal DNA of B. subtilis strain T)^\Q4:Khfq::phle [15], 
which was kindly provided by Dr. Sabine Brand. The deletion was 
accomplished by a double crossover of the antibiotic marker. B. 
subtilis strain GP1067 was constructed by transformation of the 
integrative plasmid pGP1331A/yshort into B. subtilis strain 168. 
Strain GP1067 encodes a chromosomal hfqrjlag fusion gene (C- 
terminal Plag-tag), which is controlled by the authentic promoter. 
The chromosomal insertions were verified by means of PGR and 
subsequent DNA sequencing. 

B. subtilis was transformed with DNA according to a two-step 
protocol as described [63]. Transformants were selected on agar 
containing chloramphenicol (5 Hg/ml), spectinomycin (150 Hg/ 
ml) or phleomycin (2.5 Hg/ml). 

Construction of Plasmids 

The oligonucleotides used for cloning are listed in Table S4 in 
File SI. The plasmids used in this study are based on the Gram- 
positive - E. colt shutde vector pAD123 [64], provided by the 
Bacillus Genetic Stock Center (BGSC). To remove the Shine and 
Dalgarno sequence of the ^;mut3a gene of plasmid pAD123, a 
DNA fragment was generated by means of PGR using the forward 
primer V67, which contained an Xhal site, and the reverse primer 
W67, which contained a HindLll site. The PGR product was 
cleaved with Xbal and Hindlll, and then ligated into the 
C()rr('sp()nding sites of plasmid pAD123 resulting in plasmid 
pADARBS. 

The plasmids pADresAts and pADresDts, bearing transcrip- 
tional resA-gfp and resD-gfp fusions, are derivatives of plasmid 
pAD123. Gonstruction of plasmid pADresAts: the forward primer 
F81, which contained an EcdSl site, and the reverse primer 188, 
which contained a Xbdl site, were used to amplify the region 
spanning nucleotides — 1 0 1 to +4 with regard to the transcriptional 
start of resA. The resulting PGR fragment was cleaved with EcoSA 
and Xbal, and then cloned into the corresponding sites of plasmid 
pAD123. Gonstruction of plasmid pADresDts: the forward primer 
J88, which contained an EcdSA site, and the reverse primer L88, 
which contained a Xbal site, were used to amplify the region 
spanning nucleotides — 1 00 to +4 with regard to the transcriptional 
start of resD. The resulting PGR fragment was cleaved with EcoKl 
and Xbal, and then cloned into the corresponding sites of plasmid 
pAD123. 

The plasmids pADresAti and pADresDti, bearing translational 
resA::gfp and resD::g/p fusions, are derivatives of plasmid pADARBS. 
Gonstruction of plasmid pADresAtl: the forward primer F81, 
which contained an _EcoRI site, and the reverse primer G8 1 , which 
contained a Xbal site, were used to amplify the region spanning 
nucleotides —131 to +18 with regard to the A(+l) of the start 
codon otresA. The resulting PGR fragment was cleaved with EcdRl 
and Xhal and then cloned into the corresponding sites of plasmid 
pADARBS. Gonstruction of plasmid pADresDtl: the forward 
primer J88, containing an &oRI site, and the reverse primer K88, 
containing a Xbal site, were used to amplify the region spanning 
nucleotides —124 to -1-18 with regard to the A(-l-l) of the start 
codon of resD. The resulting PGR fragment was cleaved with 
EcoRl and Xbal and then cloned into the corresponding sites of 
plasmid pADARBS. 

The plasmid pGP1331A/9short used for construction of-B. subtilis 
strain GP1067 is based on the integrative plasmid pGP1331 [65]. 
Plasmid pGP1331 was used to abut the 3'-end of hfq with the 
sequence encoding a triple Flag-tag. The hfqrjlagp) fusion gene is 
under transcription control of the authentic promoter. Gonstruc- 
tion of plasmid pGP1331A/^short: the forward primer ML152, 
containing a BamSAl site, and the reverse primer ML151, 



containing a SaK site, were used to amplify the kfq gene. The 
resulting PGR fragment was cleaved with BamHl and Sail, and 
then cloned into the corresponding sites of plasmid pGP1331. 

Growth Competition Experiment 

B. subtilis strain 168 wt and strain 168AA/§' and the strains 
GP1067, encoding the hfq::Jhg fusion gene, and 168AA/j, 
respectively, were grown concomitantiy in GS-glucose medium 
at 37°G. The strains were inoculated at an initial GFU/ml of 
IxlO', and growth was followed over 10 days by scoring the 
colony forming units (GFU). Each day serial dilutions were plated 
on LB agar plates containing appropriate antibiotics for selection. 
To verify synthesis of Hfq in strain GP1067 during the course of 
the experiment, samples for western-blot analysis were withdrawn 
at the times indicated in Fig. IB. 

Western-blot Analyses 

The cellular levels of Hfq and S2 ribosomal protein were 
determined by quantitative immunoblotting. B. subtilis strain 
GP1067, encoding the Hfq-Flag-tag fusion protein, was grown 
in GS-glucose medium at 37°G. Samples were withdrawn at 
various times (indicated in Fig. lA and B) and the cells were 
pelleted by centrifugation. The cell pellets were resuspended in 
50 |xl lysis buffer (Ix TE buffer, 20 mM EDTA, 1 mg/ml 
lysozyme) per ml cell culture and incubated at 37°G for 20 min. 
At an ODgoo of 1.0, 2 |J.l of 20% SDS were added, and the samples 
were heated at 95°G for 5 min. Then, 1 volume of 2 x protein 
sample buffer was added and the suspension was heated for 
additional 5 min at 95°G. Equal amounts of total protein were 
separated on 12% SDS-polyacrylamide gels, and then blotted to a 
nitrocellulose mc-mbrane (GE Healthcare). The blots were blocked 
with 5% dry milk in TBS buffer, and then probed with anti-Flag 
(Roth) or anti-S2 (lab stock) antibodies. The antibody-antigen 
complexes were visualized with goat anti-rabbit immunoglobulin 
alkaline-phosphatase conjugated antibody (Sigma) using NBT 
(Nitroblue-tetrazolium-chloride; Sigma) and BCIP (5-Bromo-4- 
chloro-3-indolyl phosphate toluidine salt; Sigma) in alkaline 
phosphatase-bufier (10 mM NaCl, 5 mM MgGlj, 100 mM Tris/ 
HGl, pH 9.5) as a chromogenic substrate. 

Comparative Transcriptome Analysis 

The B. subtilis strain 168 and strain 168AA^ were grown in GS- 
glucose medium at 37°G. Total RNA from logarithmically 
growing cultures (OD500 = 0.7) and cultures grown to early 
stationary phase (ODgoo = 2.0) were prepared as follows. Total 
RNA was isolated using the Trizol method (Ambion). The samples 
were treated with DNase I (DNase I, RNase-free, Roche Applied 
Science). Next, a control PGR was performed using the aprE^ 
specific oligonucleotides M70 and N70 (Table S4 in File SI) to 
confirm complete degradation of chromosomal DNA. The 
RNAseq analysis included 2 biological and 2 technical replicates. 
The RNA was fragmented to an average length of 200—300 nt by 
incubation for 2 minutes at 94°G in 40 mM Tris-acetate pH 8.2, 
100 mM potassium-acetate and 30 mM magnesium-acetate [66]. 
The samples were first cooled on ice and then purified on a 
Sephadex G50 column (GE Healthcare). The cDNA synthesis was 
carried out using the Superscript Double-Stranded cDNA 
Synthesis Kit (Invitrogen) following the manufacturer's instruc- 
tions. The cDNA was purified using phenol/GHGls. The cDNA 
of the different samples was further processed and subjected to 
next generation RNA sequencing (NGS; lUumina platform GAIIx) 
at the Vienna Biocenter Gampus Science Support Facility (http:// 
www.csfac.at). The protocol for cDNA samples processing and 



PLOS ONE I www.plosone.org 



7 



June 2014 I Volume 9 | Issue 6 | e98661 



Bacillus subtilis Transcriptome Analysis 



RNA sequencing is available at http://www.csf.ac.at/facilities/ 
ngs/protocols. 

Bioinformatic Analyses 

The obtained sequencing reads were mapped on the genome of 
Bacillus subtilis 168 (NC_000964.3) using Segemehl (version 0.9.3) 
with default parameters [67]. For each annotated mRNA from the 
NCBI database the number of mapped reads, for each growth 
phase and separately for each of the replicates, were determined 
by counting all reads with an overlap of at least 1 nt. To identify 
transcripts with dilferent abundance in B. subtilis 168 and 168Ahfq 
during logarithmic growth or at the onset of stationary phase, a 
differential gene expression analysis was performed, using the tool 
DEseq (version 1 .5), which is part of the bioconductor packages 
[68] . A p-value (adjusted for multiple testing) below 0. 1 was set as 
threshold for significance. 

Target Prediction for sRNAs 

The B. subtilis sRNAs identified and predicted in recent studies 
[50-51] were used as a source for target prediction. The RNA 
local folding opening energy (a measure for the accessibility of the 
sequence and the binding sites) was calculated for pair probabil- 
ities using the program RNAplfold [69] of the ViennaRNA 
package with the options -W 200 -L 150 [70]. The putative 
binding sites between sRNA and mRNA targets were calculated 
using RNAplex software with the options -A -a -1 30 [71]. 

Primer Extension Analyses 

The levels oihfq mRNA and 5S ribosomal RNA were determined 
by primer extension. B. subtilis wt strain 168 was grown in CS- 
glucose medium at 37°C. Aliquots for RNA purification were 
withdrawn at the ODgoo values indicated in Fig. 1 A. The total RNA 
was purified with the Trizol method (Ambion). The primer 
extension analysis was performed using 1 5 |J,g of total RNA. After 
annealing to the A/?"Specific ['^P] 5 '-end-labelled (Amersham 
Pharmacia Biotech) oligonucleotide Q65 (Table S4 in File SI) 
and addition of AMV-Reverse Transcriptase (Promega), cDNA 
synthesis was carried out for 15 min at 45°C. The oligonucleotide 
Q_65 is complementary to nt +23 to +45 with regard to the A (+1) of 
the AUG start codon of hfq mRNA. Primer extension yielded a 
product of 80 nt in length. As a loading control, the 5S rRNA levels 
were determined in the corresponding RNA samples by primer 
extension using the 5S rRNA-specific [^^P] 5 '-end-labeled oligo- 
nucleotide B64 (Table S4 in File SI). The oligonucleotide B64 is 
complementary to nt +58 to +79 with regard to the transcription 
start site (+1) of 5S rRNA. Primer extension yielded a product of 
79 nt in length. The extension products were resolved on a 8% 
polyacrylamide-8 M urea gel and were then visualized with a 
Phosphorlmager (Molecular Dynamics). 

Northern-blot Analyses 

The steady-state levels of the sRNAs 13, FsrA and SurA were 
determined by Northern-blot analyses. B. subtilis strain 168 and 
strain 168AA/g were grown in CS-glucose medium at 37°C. 
Aliquots for RNA purification were withdrawn at the ODgoo 
values indicated in Fig. 5B. RNA extraction was performed using 
the Trizol method (Ambion). 15 |.lg of total RNA were denatured 
for 5 min at 85 °C in RNA loading dye, separated on 8% 
polyacrylamide-8 M urea gels, and then transferred to a nylon 
membrane (Amersham Hybond-N) by electroblotting. The RNA 
was cross-linked to the membrane by exposure to UV light. 
Target-specific [y-'^P] -5 '-end-labeled (Amersham Pharmacia 
Biotech) ofigonucleotides were used as indicated in Table S4 in 



File SI. sRNA 13 was detected using an internally radiolabelled 
probe synthesized by means of PGR in the presence of [of-^^P]- 
dATP (Amersham Pharmacia Biotech) using primers M80/N80 
(Table S4 in File SI). The hybridization signals were visualized 
using a Phosphorlmager (Molecular Dynamics). The experiments 
were performed in dupficate. 

Gfp-assays 

B. subtilis strains 168 and 168AA/^ bearing either plasmid 
pADresAts, pADresDts, pADresAtl or pADresDtl were grown in 
GS-glucose medium at 37°G to an ODgoo of 0.7, at which time 
1 ml aliquots were withdrawn in duplicate. The samples were 
pelleted by centrifugation and the pellets were dissolved in 2 ml Ix 
TE buffer to measure the fluorescence with a spectro-fluorometer 
{Jasco FP-6300) at wavelengths Ex-501 nm and Em-511 nm. The 
experiment was performed in tripficate. 

Supporting Information 

Figure SI Growth and sampling of the B. subtilis 
strains for RNA-seq analysis. Bacillus subtilis strains 168 wt 

(triangles) and 168AA/^ (squares) were grown in CS-glucose 
medium at 37°G. Samples for RNA extraction were withdrawn 
at an ODggo of 0.7 and 2.0 as indicated by the arrows. 
(TIF) 

Figure S2 Schematic representation of the resABCDE 
operon and regulation by ResD of anaerobic respiration 
in B. subtilis. Upon autophosphorylation ResE transfers a 
phosphate to ResD. Subsequentiy, phosphorylated ResD activates 
transcription at target promoters. Activated genes include besides 
others fnr, encoding the anaerobic transcriptional regulator, 
nasDEF, which constitute an operon encoding the subunits of 
nitrite reductase', nasBd, nitrate reductase, hmp, flavohemoglobin 
and narK, nitrite extrusion protein. Shaded boxes: genes involved 
in anaerobic regulation; Arrows: regulatory flows; Arrows above 
boxes denote start sites and the direction of transcription; Red 
stars indicate an increased transcript abundance in the B. subtilis 
l68Ahfq mutant (see Tables SI and S2 in File SI). Adapted from 
Nakano and Zuber (1998) Am. Rev. Microbiol. 52: 165-190. 
(TIF) 

Figure S3 The steady state levels of type I TA system 
RNAs are reduced in the absence of Hfq. Bacillus subtilis 
strains 168 wt (wt) and l6'SAkfq (A) were grown in CS-glucose 
medium at 37°G. Samples for RNA extraction were withdrawn at 
the ODgoo values indicated on top. RNA extraction was 
performed using the Trizol method (Ambion). 15 |j.g of total 
RNA were denatured for 5 min at 85"C in RNA loading dye, 
separated on 8"/o polyacr}'lamide-8 M urea gels, and then 
transferred to a nylon membrane (Amersham Hybond-N) by 
electroblotting. The RNA was cross-linked to the membrane by 
exposure to UV light. The membrane was hybridized with target- 
specific [y-'^^P] -5 '-end-labeled (Amersham Pharmacia Biotech) 
oligonucleotides as indicated in Table S4 in File SI. The 
hybridization signals were visualized using a Phosphorlmager 
(Molecular Dynamics). Only the relevant sections of the 
autoracUographs are shown. 
(TIF) 

File SI Supplemental Tables. Transcripts displaying altered 
abundance in the absence of Hfq during logarithmic growth 
(Table SI) and in early stationary phase (Table S2). Bacillus 
subtihs strains used in this study (Table S3). Oligonucleotides used 
in this study (Table S4). 
(DOG) 
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